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£_) ' We study the finite-size and surface effects on the thermal and spatial behaviors of the magneti- 

sation of a small magnetic particle. We consider two systems: 1) A box-shaped particle of simple 
cubic structure with either periodic or free boundary conditions. This case is treated analytically 
using the isotropic model of D-component spin vectors in the limit D — > oo, including the magnetic 
field. 2) A more realistic particle (7-Fe20a) of ellipsoidal (or spherical) shape with open boundaries. 
The magnetic state in this particle is described by the anisotropic classical Dirac-Heisenberg model 
including exchange and dipolar interactions, and bulk and surface anisotropy. This case is dealt 
with by the classical Monte Carlo technique. 
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I. INTRODUCTION 



From the physical point of view, nanoparticles exhibit such interesting features as superparamagnetism and expo- 
nentially slow relaxation rates at low temperatures due to anisotropy barriers. However, the picture of a single-domain 
^- , magnetic particle where all spins are pointing into the same direction, leading to coherent relaxation processes, ceases 
to be valid for very small particles where surface effects become really crucial. For instance, in a particle of radius 

ON , ~ 4 nm, 50% of atoms lie on the surface. Therefore, it is necessary to understand the effect of free boundaries first on 
the static and then on the dynamical properties of nanoparticles. However, one of the difficulties which is inherent to 
systems of round (spherical or ellipsoidal) geometries, consists in separating surface effects due to symmetry breaking 
of the crystal field on the boundaries and the unavoidable finite-size effects caused by using systems of finite size. In 

qs hypercubic systems, this problem is easily handled by using periodic boundary conditions, but this is not possible in 

•">; . other topology, and thus surface and finite-size effects are mixed together. 

In this article, we discuss surface and finite-size effects on the thermal and spatial behaviours of the intrinsic 
magnetisation of an isolated small particle. We consider two different systems: 1) A cube of simple cubic structure 
with either periodic or free boundary conditions. This system is treated analytically by the isotropic model of D- 
component spin vectors in the limit D — ► oo, in magnetic field W. 2) The second system, which is more realistic, is 
the maghemite particle (7-Fe20a) of ellipsoidal (or spherical) shape with open boundaries. The appropriate model is 
the anisotropic classical Dirac-Heisenberg model including exchange and dipolar interactions, and taking account of 
bulk and surface anisotropy. On the contrary, this system can only be dealt with using numerical approaches such 
as the classical Monte Carlo technique M . In the case of a cubic system we obtain the thermal behaviour of local 
magnetisations at the center of faces, edges and corners. An exact and very useful relation between the intrinsic 
magnetisation and the magnetisation induced by the magnetic field, valid at all temperatures and fields, was obtained 
in Ref. (1). It was shown that the positive contribution of finite-size effects to the magnetisation is lower than the 
negative one rendered by boundary effects, thus leading to a net decrease of the magnetisation with respect to the 
bulk. For the maghemite, this study has been performed in a very small and constant magnetic field; the surface shell 
is assumed to be of constant thickness and only the particle size is varied. So, the thermal behaviour of the intrinsic 
magnetisation is obtained for different particle sizes B] . This behaviour is compared with that of a cubic maghemite 
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particle with periodic boundary conditions but without anisotropy. In this case the contributions of finite-size and 
surface effects lead to the same results as for the cube system, but the difference between them is now much larger, 
due to surface anisotropy. In addition, we show that the magnetisation profile is temperature dependent. 

II. CUBIC SYSTEM: D -> oo SPHERICAL MODEL 

We consider an isotropic box-shaped magnetic system of volume j\f = L 3 , with simple-cubic lattice structure, and 
nearest-neighbour exchange coupling, in a uniform magnetic field. For this we use the Hamiltonian of the isotropic 
classical /^-component vector model ||, that is, 

1 D 

% = ~h" / j s i ~ 2 / j^V / < s ai s aji (1) 

i i,j a— 1 

where Sj is the normalized D-component vector, |s,| = 1; h= H/Jo is the magnetic field, and Ay = Jij/Jo the 
exchange coupling. We also define the reduced temperature 9 = T/T^ FA , T^ IFA = Jq/D being the Curie temperature 
of this model in the mean-field approximation, Jo is the zero-momentum Fourier component of J y - . In this model, 
the magnetisation m is directed along the field h, so that h —he z and m^m^e^. Using the diagram technique for 
classical spin systems [pi in the limit D — > oo, generalizing it so as to include the magnetic field and adopting a matrix 
formalism jl|, one ends up with a closed system of equations for the average magnetisation component m; = (s Z i) 
and correlation functions s y = D (s a iS a j) with a > 2 H, 

J2 %™j = Gih, J2 v a s J l = 9GAl ' ( 2 ) 

3 3 

where T>ij = Sij — GiXij is the Dyson matrix of the problem, and Gi is a local function to be determined from the set 
of constraint equations on all sites i = 1 , . . . , j\f of the lattice 

su + m? = 1. (3) 

Now, we define the induced average magnetisation per site by 

m =^E m > ( 4 ) 

i 

which vanishes for finite-size systems in the absence of magnetic field due to the Golstone mode associated with global 
rotations of the magnetisation. On the other hand, it is clear that at temperatures 6 <C 1 the spins in the system are 
aligned with respect to each other and there should exist an intrinsic magnetisation. The latter is usually defined for 
finite-size systems as 



M = 




N\V*V^/"\ 



1 u 
Jp Z^i Sl i> ( 5 ) 



where the second equality is valid in the limit D — > oo. Note that M > m and that M remains non zero for h = 0; 
in this case in the limit 9 — * 0, Sy- = 1 for all i and j, and M — * 1. For 9 — * oo the spins become uncorrelated and 
M — ► 1/yW. In the limit of j\f — > oo, the intrinsic magnetisation M approaches that of the bulk system. In the 
presence of a magnetic field, the Goldstone mode is suppressed and the magnetisation m of Eq.(H) no longer vanishes, 
this is why we call it the supermagnetisation, in contrast with the intrinsic magnetisation M. If the field is strong the 
magnitude of the supermagnetisation approaches the intrinsic magnetisation. 
An important exact relation was established in Ref. jl| between M and m, 

= M 2jVMh/9 = MB{NMH/T) (6) 

1 + y/1 + {2NMh/9f 

where B(£) = (2£/D)/ 1 + y/l + (2£,/D) 2 is the Langevin function for D > 1. Note that Eq. @ is usually applied 
to superparamagnetic systems with the spontaneous bulk magnetisation m\,(T) in place of M{T 1 H). However, unlike 



rrih(T), the intrinsic magnetisation M of Eq. (g) is a pertinent characteristic of a finite magnetic system and depends 
on both field and temperature. 

Solving the model above consists in determining m^ and Sij as functions of Gi from the linear equations (En, and 
inserting these solutions in the constraint equation (j3|) in order to obtain Gi . Two types of boundary conditions are 
considered, free boundary conditions (fbc) and periodic boundary conditions (pbc). In the case of fbc, m^ and Gi are 
inhomogeneous and Sy non-trivially depends on both indices due to boundary effects. In this case the exact solution 
is found numerically, though some analytic calculations can be performed at low temperature and field. Whereas in 
the pbc case the solution becomes homogeneous and the problem greatly simplifies. Although the model with pbc is 
unphysical, it allows for an analytical treatment and study of finite-size effects separately from boundary effects. 

At low temperature, the intrinsic magnetisation in the fbc case, including only the contributions from faces, reads 
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where W is the well-known Watson's integral and 



Aat = 



^r^S 1 - 1 ^] 



>0 



(8) 



describes the finite-size effects, with Ajv ex 1/L, while the last term in (pf) represents the contribution from boundaries. 
The first term, on the other hand, is the bulk contribution which survives in the limit L — > oo. In contrast with the 
finite-size effects, boundary effects entail a decrease of the intrinsic magnetisation. The contributions to Eq. (M) from 
the edges and corner are of order 9/L 2 and 0/L 3 , respectively. 
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FIG. 1. Temperature dependence of the intrinsic magnetisation M, Eq. (0), and local magnetisations of the I4 3 cubic 
system with free and periodic boundary conditions in zero field. 

Fig. yj shows the temperature dependence of the intrinsic magnetisation M , Eq. (|3|), and local magnetisations of 
the I4 3 cubic system with free and periodic boundary conditions in zero field. For periodic boundary conditions, M 
exceeds the bulk magnetisation at all temperatures. In particular, at low temperatures this agrees with the positive 
sign of the finite-size correction to the magnetisation, Eq. ([?]). The magnetisation at the center of the cube with free 
boundary conditions is rather close to that for the model with pbc in the whole temperature range and converges 
with the latter at low temperatures. Local magnetisations at the center of the faces and edges and those at the 
corners decrease with temperature much faster than the magnetisation at the center. This is also true for the intrinsic 
magnetisation M which is the average of the local magnetisation Mi over the volume of the system. One can see 



that, in the temperature range below the bulk critical temperature, M is smaller than the bulk magnetisation. This 
means that the boundary effects suppressing M are stronger than the finite-size effects which lead to the increase of 
the latter, and this is in agreement with the low-temperature formula of Eq. ([]). 

III. MAGHEMITE PARTICLES: MONTE CARLO SIMULATIONS 

In this section, we consider the more realistic case of (ferrimagnetic) maghemite nanoparticles (7-FC2O3) of ellip- 
soidal (or spherical) shape with open boundaries, in a very small and uniform magnetic field. The surface shell is 
assumed to be of constant thickness (~ 0.35 nm) in J4|, and only the particle size is varied J2j. 

To deal with spatial magnetisation distributions |5| one has to consider exchange, anisotropy and magneto-static 
energies together. Accordingly, our model for a nanoparticle is the classical Dirac-Heisenberg Hamiltonian including 
exchange and dipole-dipole interactions, anisotropy, and Zeeman contributions. Denoting (without writing explicitly) 
the dipole-dipole interaction by Hdip, our model reads, 

_ N t N t 

n = - E J a ? E s ? ' s f+n - ^E ( s < • e <) 2 - ^b)h E s, + H dip 

i,i\ a.J3 i= 1 i— 1 

where J a p are the exchange couplings between nearest neighbours spanned by the unit vector n; Sf is the (classical) 
spin vector of the a th atom at site i; H is the uniform field applied to all spins in the particle, K > is the anisotropy 
constant and e^ the single-site anisotropy axis. In both cases of a spherical and ellipsoidal particle, we consider a 
uniaxial anisotropy in the core along our z reference axis (major axis for the ellipsoid), and single-site anisotropy on 
the surface, with equal anisotropy constant K s , and ej are defined so as to point outward and normal to the surface 

Our method of simulation proceeds as follows: we start with a regular box of spinel structure, then we cut in 
a sphere or an ellipsoid that contains the total number Nt of spins of a given particle. We distinguish between 
spins in the core (of number N c ) from those on the surface (N s ) of the particle according to whether or not their 
coordination number is equal to that of a system with periodic boundary conditions (pbc). All spins in the core 
and on the surface are identical but interact via different couplings; exchange interactions between the core and 
surface spins are taken equal to those inside the core. Our parameters are as follows: the exchange interactions 
are (in units of K) JabI^-b — —28.1, Jbb/^b — —8.6, JaaI^-b — —21.0. The bulk and surface anisotropies are 
k c = {Kc/ks) — 8.13 x 10~ 3 , k s = (Ks/ks) — 0.5, respectively, where ks is the Boltzmann constant. 

In Fig. @, we plot the thermal variation of the core and surface contributions to the magnetisation (per site) as a 
function of the reduced temperature r core = T/T™ re for N t = 909,3766 corresponding to N st = N s /N t = 53%, 41% 
and diameters of circa 4 and 6 nm, respectively. The core and surface magnetisations are averages over all spins in 
the core or on the surface, respectively. For both sizes we see that the surface magnetisation M sur f decreases more 
rapidly than the core contribution M c as the temperature increases, and has a positive curvature while that of M c 
is negative. Moreover, it is seen that even the (normalised) core magnetisation per site does not reach its saturation 
value of 1 at very low temperatures, which shows that the magnetic order in the core is disturbed by the surface (see 
Fig. 4 below). As the size decreases the maximum value of M sur f decreases showing that the magnetic disorder is 
enhanced. 

In Fig. H we plot the core and surface magnetisations (with N t = 909, 3766, and N st — 53%, 41%), the magnetisation 
of a cube with spinel crystalline structure and pbc, and the bulk magnetisation as functions of the reduced temperature 
T core . Apart from the obvious shift to lower temperatures of the critical region due to the finite-size and surface effects, 
we see that, as was also shown analytically for the cube system, the finite-size effects give a positive contribution 
to the magnetisation with respect to the bulk, whereas the surface effects yield a negative contribution. Moreover, 
it is seen that for nanoparticles the contribution from the surface is much larger than that coming from finite-size 
effects. The difference between the two contributions appears to be enhanced by the surface anisotropy in the case of 
nanoparticles. 

In Fig. [I we plot the spatial evolution of the orientation of the magnetic moment from the center to the border 
of the particle, at different temperatures. At all temperatures, the magnetisation decreases with increasing particle 
radius. This obviously suggests that the magnetic disorder starts at the surface and gradually propagates into the core 
of the particle. At high temperatures, the local magnetisation exhibits a jump of temperature-dependent height, and 
continues to decrease. This indicates that there is a temperature-dependent radius, smaller than the particle radius, 
within which the magntisation assumes relatively high values. This result agrees with that of Ref. J7] (for spherical 
nanoparticles with simple cubic structure) where this radius was called the magnetic radius. The local magnetisation 
also depends on the direction of the radius vector, especially in an ellipsoidal particle. 
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FIG. 2. Temperature dependence of the surface and core magnetisations for N t = 909 and 3766. 




FIG. 3. Temperature dependence of the surface and core magnetisations for Nt — 3766, magnetisation of the bulk system 
(Nt — oo), and the magnetisation of the cube with the spinel structure and periodic boundary conditions (PBC) with Nt = 40 3 . 




Particle radius 

FIG. 4. Spatial variation of the net magnetisation of a spherical nanoparticle of 3140 spins, as a function of the normalised 
particle radius, for r core < 1, and r core = 0.5, r core ~ 1". 

IV. CONCLUSION 

Both for the cube system and the nanoparticle of the maghemite type surface effects yield a negative contribution 
to the intrinsic magnetisation, which is larger than the positive contribution of finite-size effects, and this results in 
a net decrease of the magnetisation with respect to that of the bulk system. In the first case we have been able 
to separate finite-size effects from surface effects by considering the same system with periodic and free boundary 
conditions. On the other hand, the results for a spherical or ellipsoidal nanoparticle with free boundaries have been 
compared to those of a cube with a spinel structure and periodic boundary conditions, but without any anisotropy. 
In this case, it turns out that the contributions from surface and finite-size effects have the same sign as before but 
the difference between them becomes larger, due to surface anisotropy. 

These spin models invariably predict that the surface magnetisation (per spin) of systems with free boundaries 
is smaller than the magnetisation of the bulk system. However, experiments on layered systems, especially of 3d 
elements, have shown that there is enhancement of the magnetic moment on the surface, which has been attributed 
to the contribution of orbital moments ||. It is clear that the models presented here do not account for such effect, 
but they can be generalised so as to include orbital as well as spin vectors. 
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